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Abstract 

We compute fourth sound for superfluids dual to a charged scalar and a gauge field 
in an AdS4 background. For holographic superfluids with condensates that have a large 
scaling dimension (greater than approximately two), we find that fourth sound approaches 
first sound at low temperatures. For condensates that a have a small scaling dimension 
it exhibits non-conformal behavior at low temperatures which may be tied to the non- 
conformal behavior of the order parameter of the superfluid. We show that by introducing 
an appropriate scalar potential, conformal invariance can be enforced at low temperatures. 
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1 Introduction 

When '^He is cooled down to temperatures below Tx ~ 2. 17 ° K it undergoes a phase transition 
from a liquid phase, helium 1, to a superfluid phase, helium 11. While helium 1 behaves like 
a normal low-viscous fluid, helium 11 exhibits both viscous and non-viscous behavior. This 
property of helium 11 (and several other properties) can be explained by a two-fluid model 
flrst suggested by Tisza [I] and expanded upon by Landau [2j. In this model, helium 11 is 
described as a fluid with two components. One, the normal component, has normal viscous 
behavior, while the other, the superfluid component, behaves like a frictionless fluid. 

Since the gauge-gravity duality O HI [5] has been successful in describing the hydrody- 
namic regime of gauge theories with holographic duals [01 E] , it is interesting to ask to what 
extent it can describe the hydrodynamic behavior of superfluids. To create a "holographic 
superfluid," or a gravity dual of a system with a condensed phase, one needs a mechanism 
by which a condensate will be generated in an asymptotically AdS geometry. In [8] it was 
argued that a charged scalar coupled to an abelian gauge fleld in AdS space (an AdS version 
of the abelian Higgs model) may condense in the presence of a black hole horizon. From the 
point of view of the boundary theory this corresponds to a second order phase transition 
whose order parameter is dual to the condensed scalar [9] . Other types of symmetry breaking 
mechanisms were considered in [T0| [TTl [T2| [T3] . 

In [T^ [T5] some similarities between the phase diagram of the condensate dual to the 
abelian Higgs model in the bulk, and that of a superfluid were discussed. Like a superfluid, 
the phase diagram of the boundary theory was shown to include a tricritical point and a 
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critical superfluid velocity. A similar analysis was carried out in [131 IIS] fo^' ^ holographic 
p-wave superfluid derived from a non-abelian gauge field in the bulk. 

In this work we investigate the phase velocity of particular sound modes in holographic 
superfiuids. In one component fluids, small perturbations of the pressure follow a wave 
equation with a phase velocity equal to the speed of sound. In superfluids, there are two 
propagating modes involving the pressure and the entropy density. As a consequence, su- 
perfluids support two sound modes and, appropriately, two phase velocities called flrst and 
second sound. In the limit where the density of the superfluid component goes to zero, sec- 
ond sound vanishes and the remaining sound mode is flrst sound — the one associated with 
pressure waves. Another limit can be obtained when an external force prevents the normal 
component from flowing. In this case a propagating mode involving density fluctuations will 
still remain. Since the normal component is stationary, these density fluctuations necessarily 
involve non-vanishing superfluid velocity. The phase velocity of such modes is called fourth 
sound [I7]E1 See also [TS]. For ^He its value coincides with second sound near the phase 
transition, and flrst sound at low temperatures [17] . 

Going back to a holographic setup, in [15] the authors found an approximate expres- 
sion for second sound, valid at high temperatures and the probe approximation (where the 
geometry does not back-react on the matter flelds). In the following section we show that 
one may reinterpret the approximate expression obtained in [15] as an exact expression for 
fourth sound, and then extend its range of validity to low temperatures, as long as the probe 
approximation holds. 

From the results of [T^] one flnds that, in the abelian Higgs model, at low temperatures, 
fourth sound squared approaches 1/2 for dimension two condensates and 1/3 for dimension 
one condensates in d = 3 space-time dimensions. In section [2] we explain why, in a conformal 
theory, fourth sound should approach flrst sound at low temperatures 



where d is the spacetime dimension. Equation ([T]) explains the asymptotic value of fourth 
sound found in [15] for dimension two operators but presents a puzzle for dimension one 
operators. In sections [3] and H] we use the Abelian Higgs model of [8],|9] to study holographic 
superfluids with condensates of dimension 1/2 < A < 5. We flnd that for 1/2 < A ^ 2 
conformal symmetry is apparently broken at low temperatures; the breaking of conformal 

^ Third sound, also defined in [17j , is the phase velocity of wave propagation on the surface of a thin film 
of superfluid. 
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symmetry introduces an extra scale which allows the asymptotic value of fourth sound to 
deviate from its expected conformal value. As we show in section H] this extra scale can 
be tied to the divergence of the order parameter at low temperatures. Such a divergence, 
for dimension one operators, was first discussed in [9]. In section [5] we adjust the matter 
action so that the scalar potential is bounded from below, enforcing conformal behavior at 
low temperatures. We summarize our findings in section O 

2 Sound modes in the Tisza-Landau two-fluid model 

Hydrodynamics can be thought of as an effective theory, valid as long as the relevant fields 
vary slowly in space and time relative to some microscopic scale. In the absence of sponta- 
neously broken symmetries, the hydrodynamic variables are the temperature, T, chemical 
potential, /i, and fluid velocity, u'" with m = 0, . . . ,d — 1. In the presence of a sponta- 
neously broken symmetry, the Goldstone boson provides an extra degree of freedom. The 
hydrodynamic variables in such a system are T, /i, and ^"^ = d"^ip. In [12] it was shown 
that, at the inviscid level, one may consistently construct a conserved energy momentum 
tensor and a conserved current out of these variables. Further, the resulting formulation of 
fluid dynamics coincides with the relativistic version of the Tisza-Landau two-fluid model of 
[2D1 EH- See also appendix A of [B]. 

In the formulation of [19] the energy momentum tensor, T"*", takes the form 

j.mn = (e + p) y'"^" + Prj^^n + /^Cr (2) 

where the energy density, e, and the pressure, P, are related to the entropy, s, temperature, 
T, density of the normal phase, pn, and chemical potential, /i, through 

e + P = Ts + /ipn . (3) 

The density of the superfluid component is denoted ps = and the superfluid velocity is 
given by f ™ = / p. For the current one has 

J'" = pnM™ + PC" ■ (4) 
In addition, the superfluid velocity satisfies a "Josephson equation" 

U'^im = -/^ (5) 
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and we remember that it is generated by the gradient of the Goldstone boson y?, 

Cm = . (6) 

In the absence of sources, the equations of motion are given by 

-g^T'"" = (7) 
9„J"^ = 0. (8) 

It is often useful to decompose ([7]) into a component orthogonal to the velocity field, and a 
component parallel to it. For the component parallel to the velocity field, we find from ([6]), 
(ED and 

dP = sdT + p^dfi - fCdU (9) 

that 

= UrndnT''^ = d^{su^) . (10) 

This implies that the superfluid component does not carry entropy. 

To study sound modes of the superfluid, we look at linear perturbations of the hydrody- 
namic variables around a static configuration of the fluid. In what follows, we will denote 
the unperturbed configuration with a subscript and the perturbed configuration with a 5. 
Thus, = 1, = 6v\ vP = 1, = ps = Pso + '^Ps etc. The index i runs over the spatial 
coordinates i = l,...,(i— 1. At the linear level, ©, ([6]), ([7]) and ([8]) take the form 



dt5s + sodi5u' = Q (11a) 

(eo + Po)dM + ^^oPsA^v' = -d,6P (lib) 

dt{Sps + 5pn) + pnodiSu' + psodi6v' = (11c) 

fiodM = di6p. (lid) 



The first two equations correspond to the linearized version of (fTOj) and the remaining compo- 
nents of ([7]), the third equation corresponds to the linearized current conservation equation 
([H]) and the last equation follows from and (jSD • It is possible to recast (ITT]) as wave 
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equations for the total density p = Ps + Pn and the entropy per particle a = s/p, 

dfSa = ( d.d^5P - ^"-W g^g^^^") (12) 



Po wq V Po 

Q^Sp =P^( Pn^d.d^SP + ^^^P^d,d'6p) (13) 
Wo V Po Po / 

where we have denoted Wq = eo + Pq. If we now express the chemical potential and pressure 
in terms of the entropy per particle and density, we get a set of coupled wave equations. 

Consider the case where the superfiuid phase is absent, pso = 0- Then, the entropy does 
not oscillate and we get a density wave whose phase velocity squared is 

• (14) 

This is the expected speed of sound in normal fluids. Once ps two coupled entropy- 
density waves are allowed. The phase velocities of the corresponding modes are the two 
roots of a second order polynomial. The bigger root reduces to (fT4|) when ps is set to zero. 
The smaller root is the phase velocity of a new mode called second sound. To measure second 
sound one can generate temperature waves in a tube filled with superfiuid helium, and check 
when the wavelength of the mode matches the system size and generates a resonance. The 
first measurements of second sound in helium II are described in [22] . 

Now consider a setup where one allows only the superfiuid component to flow. Such a 
situation can be achieved by scattering the normal, viscous, phase off stationary scatterers. 
In this case momentum is not conserved due to the force required to keep the scatterers 
in place. In such a setup we would set 6u^ = and remove the momentum conservation 
equation (lllbp from (fTTI) . Instead of (fT2|) and (fT3l) we would have 

d'^Ss = (15) 

d^6p = f^did'p. (16) 
Po 

Thus, we have one sound mode, similar to first sound since it excites a density wave, but 
with vanishing normal-component velocity. The existence of such a sound mode was first 
realized in [17] and is called fourth sound. From (fT6|) we find that the phase velocity of 
fourth sound is given by 

vl = ■ (17) 
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Measurements of fourth sound can be carried out by allowing the superfluid to flow through 
a very narrow tube packed with fine powder which prevents the normal component from 
flowing. This is called a superleak, since it is only the super-component of the fluid that 
leaks into the tube. Transducers are then placed at the ends of the tube to receive and 
source the sound waves [23] . 

In [T3] , equation f[T7l) was used as an approximation to second sound at high temperature 
(or large entropy). At low temperatures most of the fluid will be in the superfluid phase so 
we can approximate ps ~ p. Also, in a conformal fluid, where the only remaining scale at 
T = is /i, we expect that p oc p^^^ . Thus, for conformal fluids at low temperatures, fourth 
sound coincides with normal sound, 

vi — > . m 

T-*o d-1 

This behavior is similar to that of fourth sound in ^He, where, as mentioned earlier, it interpo- 
lates between second sound close to the phase transition and first sound at low temperatures 
[ni[23]. 



3 Setup 

Our construction of a holographic superfluid is very similar to that of [T3l [T5] . The bulk 
theory is the Abelian Higgs model studied in [SI [HI [211 [2S1 [21] , whose action is given by 

S = J (CeH - -^^mattcr) V^d^X (18) 

with 

Ceh = R-A, ^matter = \d^^ " ^qA^^\^ + ^(1^^ + ^F^.F^^ (19) 

where F^j^ = df^A^, — d^A^ and our units are such that A = —6. We will work in the probe 
approximation where the charge of the scalar field is very large, and the matter content of the 
theory decouples from gravity: Setting \& — > — > v4^/g, and scaling the coefficients 
in l^d^n such that V{\^'^) V{\^'^)/q^, implies ^^atter ^ A^attcr/g^- In the g ^ oo 
limit £mattcr dccouples from £eh- As noted in [TI], for large but finite q this approximation 
is valid only as long as \1/ and $ are not too large. In principle, for any T > 0, one may 
choose a large enough q so that the approximation will be good enough. Another way to 
think of the probe approximation is to consider a formal expansion of the full back-reacted 
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geometry in inverse powers of q. The leading order matter solutions are while the 

leading order metric is 0{q^) and receives corrections. 

Focusing first on C-^ni the solution to Einstein gravity in 3 + 1 dimensions in the presence 
of a negative cosmological constant is given by the AdS-Schwarzschild black hole, 

g^^dxV = ^ (^-f{z)dt^ + dx^, + , (20) 



where 



m = 1-3 (21) 

^0 



and /i = 0, . . . , 3. The asymptotically AdS boundary is located at ^ = and the horizon is 
located at z = zq. The temperature of the black hole, T, is given by zq = S/AttT. Using 
\E' = with ip real, Aq = (p and Ai = ai (sometimes denoted a), the equations of motion 
for the matter fields in the background fl20l) in the probe approximation are 



Sl^-^ (22b) 
S.(/8a)=^. (22c) 

Z^ 

These configurations have been studied in [H [9l [HI [271 [IS] for 

y(l^l) = m^l^l^ (23) 

and = —9/4, —2, and 0. Using the standard AdS/CFT mapping, [U [5] the mass 
of the scalar field, ip, is related to the dimension A of its dual operator (Oa), through 
= A(3 - A). Unitarity implies that A > 1/2. 

The boundary conditions we apply to f[22|) are that ip and a are finite at the horizon and 
that vanishes there. This fixes three of the six integration constants in the solution to 
f l22p . Near the asymptotically AdS boundary we require that 0(0) and a(0) are finite and 
that the z^~^ term in a near boundary series expansion of ip vanishes. 

According to the AdS/CFT dictionary [11[5],[28], the coefficient of the z^ term in a near 
boundary series expansion of ip is proportional to (Oa) while the coefficient of the z^~^ 
term is the source associated with (Oa). Our boundary conditions were chosen so that the 
source vanishes. The chemical potential fi and the charge density p of the dual theory can 
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be read off a near boundary expansion of 0, = jj, — pz + . . .. Similarly, in the gauge we are 
using, the velocity field ^ and the spatial part of the current J are given by the expansion 
a = — Jz + . . .. See [I5] for details. 



4 Solving the equations of motion 

One solution to fl22|) - fl23|l with the boundary conditions discussed in the previous section is 



/i 1- 



= 0, 



0, 



(24) 



Since (Oa) = 0, this solution describes a non-condensed phase. From the analysis of [8] we 
expect that at low enough temperatures other solutions to fl22l) - fl23p will exist, where the 
scalar field, ip, is non vanishing. Such solutions usually differ by the number of nodes of ip. 
Since solutions with nodes are expected to be unstable we will focus on the solution where 
%l) has no nodes. This will describe the condensed phase. At some critical temperature, T^, 
the solution corresponding to the non-condensed phase and the solution corresponding to 
the condensed phase overlap: There exists a non trivial ip which solves f l22al) with and a as 
in (l24l) . which is determined only up to a multiplicative integration constant — a zero models 
Computing numerically is straightforward. It involves searching for a zero mode of the 
linear equation fl22al) with = p{l — z/zq) and a = 0. In figured] we plot the dependence 
of the critical temperature on A as obtained from such a numerical analysis. A similar plot, 
which goes beyond the probe approximation we are considering here, can be found in |29j . 
Thus, at least for T < two solutions to fl221) exist: a non-condensed solution given by ([23]) 
and a solution with a non vanishing expectation value for the condensate corresponding to 
the superfluid phase. If the transition is first order, the superfluid phase may also exist at 
T>T,. 

Following [in], to see whether it is the superfluid phase or the non-condensed phase which 
is preferred, we need to compute the free energy Q given hj Q = — TS'totai where T is the 
temperature and S'totai is the total renormalized matter action. This renormalized action is 
composed of two terms. The first is the on-shell matter action f]T8]) . 



5, 



V 



on— shell 



lim (-0(9^0 + fd ■ d^d + fi)d^%l)) 



20 



^2(_02/+|^|2) 



dz 



(25) 



^ In the case of symmetry breaking in a non-abelian setup in AdSs , an analytic expression for the zero 
mode was found in [13! • I'^ |16j it was shown that information about the phase transition can be obtained 
from this zero mode without the need to resort to numerics. 
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where V is the volume of the transverse dimensions. The second contribution to Stotai 
comes from holographically renormahzing the theory [301 El], i-e., from adding boundary 
counterterms, ^counter, to ^on-shcu which render f l25|) finite. After inserting a near boundary 
expansion of a^, and ip into fl25l) . we find that in the presence of no n- vanishing sources 
for (Oa), the only divergent contributions to ^on-shcU come from the z~'^ipdzip term on the 
right hand side of (l25l) . Taking into account that the source terms for irrelevant operators 
should be treated as infinitesimal [5l [32] , we find that only a mass counterterm is required 
to keep S'on-shcii finite (meaning, a counterterm proportional to ip"^). Since this counterterm 
contribution vanishes for the boundary conditions we are interested in, we are free to ignore 
it when computing Q. For all values of A we consider in this work, we found that the 
superfluid phase is preferred over the normal phase, and that the free energy is a smooth 
function of the temperature indicating a second order phase transition. 

We used a shooting algorithm to compute the solution to fl22|l -(!23 |) corresponding to the 
superfluid phase. Instead of shooting from the boundary, we fixed dz4>{zQ) and 0(^0) and 
searched for a value of ip{zo) for which the near boundary z^~^ coefficient in a series expansion 
of ip vanishes. We did this once for 0(2:0) = and then once more with a small enough value 
of a(2o) to obtain ps = /ilim^^o l^l/l^l fo a good approximation. This procedure was carried 
out for A between 13/25 and 5, with steps ranging from 1/10 for A close to a half, to 1/2 at 
larger A. The resulting values of the speed of fourth sound, computed from ffTTl) . for A < 2 
are shown in figure [2j For A > 2, the curve for fourth sound closely follows that of A = 2. 
From figure [2] we see that for A ^ 2 and the range of temperatures we are considering, the 
asymptotic values of f| deviate from the expected conformal behavior ([T]). The reason for 
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Figure 2: (Color online) Fourth sound versus temperature for condensates of dimension 
2/3 < A < 2 in the abelian Higgs model. The temperature is measured relative to the 
critical temperature T^.. The curves are color coded according to the values of A. 
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Figure 3: Values of oa and tt-a, obtained by fitting the total charge density p = oa/U^ (f)"^ 
to the numerics. The fit was carried out for data points for which [p — ps) / p < 10^^ where 
ps is the charge density of the superfluid phase. A is the dimension of the condensate. 



this is the non conformal behavior of p at low temperatures. Numerically, we find that for 
(p — Ps)/p < 10~^ one can approximate 

P = aA/i'(|)"" (26) 

where a a and ua are plotted in figure [31 Since riA > for A ^ 2, the T ^ equation of 
state relating fi and p deviates from the expected conformal one. 

In some instances, it is difficult to decide from figure [2] whether f| will reach its conformal 
value at T = or not; we find that for 3/2 < A < 2, the curve for vf seems to level off at 
a slower rate than that of other condensates. An extrapolation of the curve for A = 3/2 
down to T = gives vf —>■ 0.48. Also, a least squares fit of the low temperature data to (126]) 
gives a much better result than fitting, say, a power law correction to a p oc behavior. On 
the other hand, the curves for condensates of dimension 3/2 < A < 2 seem to posses an 
inflection point near O.lTc. 

In figure H] we have plotted the dependence of the ratio {O/^Y^^iJi/p on T/Tc. Since 
{OaY^^h/p approaches a constant as T ^ 0, this together with (126!) implies that {0^)/ p,^ 
diverges at low temperatures whenever A ^ 2. Thus, up to the numerical uncertainties 
discussed in the previous paragraph, the non conformal behavior of (Oa) is tied to the 
non-conformal behavior of p. 

To better understand the low temperature limit of (Oa) and its relation to p, consider 
f l^^l) with a = 0. Instead of the z coordinate we will work with a rescaled radial coordinate 
C, = zji where I is some length scale to be determined shortly. Defining also the rescaled 
variables Co = -^o/-^, O = (Oa)^^, R = pi"^ and M = p,i and rescaling = (p/i, ip = ip( we 
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Figure 4: (Color online) The dimensionless ratio {O/^Y/^ji/p versus temperature. At low 
temperatures this ratio approaches a constant, which implies via fl26l) that (Oa)//^^ diverges 
at low temperatures for small A. The curves are color coded according to the value of A. 



can formally expand and ip in large Co- To leading order, we find 

dli, = (27) 

al, . (i^ili^ - (28) 

As long as (fP"il) ^ 1 we may approximate 

Vi(C) = OC^-' (29) 
in which case the solution to ( 127|) takes the form 



2a; r(i) 51 V A 

where K is a modified Bessel function of the second kind. In [5] it was observed that for 
A = 1, expression (!29l) for ip fits the low temperature numerical solution. The boundary 
conditions for (!30|) are that 0(0) = M and that does not increase exponentially in the deep 
interior. To fully justify the latter choice of boundary conditions we would need to match 
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Approximation 



1 2 3 4 5 



Figure 5: (Color online) The asymptotic value of fi/ p at low temperatures, as a 

function of A, the dimension of the scalar condensate. The blue curve shows the numerical 
result, displayed in figure HI and the dashed black curve shows the analytic approximation 
dSI]), vahd for large (Oa)//^^- 



the large C, expansion of the solution to ( 1271) with a near horizon expansion of the solution 
to fl22b|) . Such an analysis was carried out in a somewhat different context in [33]. From 
f l30|) we can compute 

As stated earlier, (l29l) is valid only as long as the 0^ term in (129|) is negligible. From 
f l30|) this implies that M^O ^ 1. Thus, we need to choose a scale £ for which both zq/ 1 
and p'^{0/s)C~^''^ <^ 1, and which would be compatible with our low temperature analysis, 
Zofi ^1. If {0^)/ii^ diverges at low temperatures (large zq^), we can use = (Oa). 
In figure [5] we plot our approximation (I3T!) versus the numerical value of (OA)^^^/i/p for 
1/2 < A < 5. 



5 Bounding the scalar potential 

It is tempting to think of fl27|) and fl251) as the zero temperature description of the condensate. 
Recall however, that we're working in the probe approximation where the matter fields are 
small enough so that they don't distort the metric. In practice, once ip and become too 
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large the geometry will back-react on the matter fields. This will certainly have a significant 
effect whenever A ^ 2, as indicated by the divergent behavior of (Oa)//^^, or by the non 
conformal behavior of fourth sound at low temperatures. One way to ensure a conformally 
invariant state at low temperatures, valid in the probe approximation, is to consider a scalar 
potential, V, which has a minimum at a finite value of |\&|. This is exactly the situation 
considered in [21], where the equations of motion of the full gravity- matter action of (ITSl) 
were studied for a scalar potential of the form 

Vm) = m'\^P\' + ^u'\^^\'' (32) 

with = —2. In the presence of the potential (132|) there exists an empty AdS solution to 
the equations of motion whenever ip is at an extremum of V, and and a.j vanish. As shown 
in [23], this allows for a condensed phase where ip fiows from ip = near the boundary to 
%jj = in the infrared. Since the geometry is AdS4 both in the ultraviolet and in the 

deep interior, one can imagine a limit where it is approximately AdS4 throughout, ensuring 
the validity of the probe approximation near T = 0. In this case we expect that close to 
T = 0, (Oa) oc /i^ and f| 1/2. This is indeed the case, as we show below. 

In the probe approximation, the equations of motion for the gauge field and scalar field 
with a scalar potential as in (132!) can be read off of (!22l) . The boundary conditions we 
used and the numerical procedure for solving the equations of motion are similar to those 
discussed in sections |3] and HI We obtained numerical solutions for A = 1 and A = 2 and 
10~^ < < 1/3 in steps as small as 10~^. The resulting value of fourth sound for dimension 
one operators is plotted in figure [61 We find that, as expected, once u is non vanishing, fourth 
sound tends to reach its conformal value, = 1/2 at low temperatures. At intermediate 
temperatures, and for relatively small u, fourth sound seems to pass a transition region where 
f| is approximately 1/3 as was the case for u = 0. The fact that vf tends to 1/2 at low 
temperatures indicates that conformal invariance is retained. Indeed, for 2/15 < < 1/3 
we find that 

p = au-^'f/ (33) 

with a = 0.69 and n = 1.07 to a fairly good approximation. For dimension two operators 
and the same range of u, we obtain (l33l) with a = 0.28 and n = 3 x 10""^. 

Computing (Oa) for A = 1, 2 we find that, unsurprisingly, once u ^ 0, {Ois^y/^ii/ p 
asymptotes to a finite value at small T. An analytic expression for the low temperatures 
dependence of (Oa) on m, p, ^ and J can be obtained as follows. Consider (12^ with a 
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Figure 6: (Color online) Fourth sound versus temperature for a dimension one condensate 
and a potential as in fl32l) . The curves are color coded according to the values of u. As u 
increases, lower temperatures become numerically accessible. 
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potential V as in (132!) . Defining lip = ipz, and setting z^fi ^ 1, the low temperature equations 
of motion take the form 

dli) = -(jy'^iP + lap^/i + u'^^^ (34a) 
Since 2; does not appear explicitly in (!34|) . we can construct the conserved quantity, 

n = \ [{dM + \d.a? - {OA?) + \[^'- l«r - luH") . (35) 

Deep in the infrared, at very large z, we expect the scalar field to be at the minimum of 
the potential (1521) . = v^/n, and the other fields to vanish, a = and = 0. With this 
information, evaluating fl35|) at large z we find 

7i = 0. (36) 

Using (|36|1 and the near boundary expansion of and a in fl35|l . we obtain 

W = /-|J? (37) 

for dimension two operators, and 



= l^if^' - leT) \&) - - |j|2) j (38) 

for dimension one operators. If the superfluid velocity vanishes, = J = 0, we find that in 
the u limit (l37j) reduces to 

{O2) , 



P 

while the '-' branch of (l38l) reduces to 



(39) 



= 1 (40) 



which coincides with (|3TD for A = 1. We've checked that (|37j), (|38|), ([39]) and (jlO]) agree with 
our numerical results and with those in the literature [T5l. 
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6 Summary 



The results of our analysis can be summarized by figures [21 [3l H] and (6) In figure [2] we see 
that in (i = 2 + 1 dimensions fourth sound deviates from its expected conformal behavior 
at low temperatures for condensates of dimension A ^ 2, due to an anomalous scaling of 
the charge density (figure [3]) • Figure H] together with figure El relate the anomalous scaling 
of the charge density to the divergent behavior of (Oa)//^^ at low temperatures. Finally, 
in figure [6] we see how one can enforce conformal behavior at low temperatures by ensuring 
an asymptotically AdS solution in the deep interior. The relation we've found between the 
charge density and the value of the condensate seems quite robust: one may summarize it by 
saying that if {0/\) / is finite at low temperatures then fourth sound will asymptote to its 
conformal value. This behavior is in agreement with the 3 + 1 dimensional p-wave superfluid 
studied in [TB] . 

Apart from the asymptotic behavior of (Oa), described at the end of sections H] and 
our results were numerical. As such, we were unable to obtain a solution at arbitrarily low 
temperatures. For n = and condensates of dimension A ^ 2 it is difficult to study the 
low temperature solutions due to the divergent behavior of p and (Oa). For large values of 
A, keeping track of (at least) the first 2A — 3 terms in a near boundary series expansion 
of ip becomes challenging as Zq^l becomes large. Unfortunately, this numerical issue makes 
it hard to find an exact critical dimension, Ac, below which the the scaling dimension of 
the condensate becomes anomalous at low temperatures. To obtain a better estimate of A,,, 
we need to go down to temperatures of at least half the current value obtained. With our 
current algorithm, probing such low temperatures would require significantly more precision 
than the 40 digits of precision we have worked with. 

Experimental results for fourth sound in ^He [23] closely follow the theoretical curve 
predicted in [T7j. It is satisfying that for condensates with A ^ 2 the dependance of fourth 
sound of holographic superfluids on the temperature is very similar to that of helium II. With 
such a qualitatively good fit to fourth sound, it would be interesting to compute second sound 
by solving the equations of motion for the full back-reacted geometry. 
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